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Axon guidance involves multiple second messenger signal transduction pathways. Although each signal 
transduction pathway has been characterized, only a few studies have examined crosstalk between these 
cascades. Here, we applied a simultaneous second messenger imaging method to the growth cone and 
demonstrated correlations between cAMP, cGMP, and Ca 2+ . The levels of cAMP and cGMP in 
non-stimulated freely extending growth cones showed a negative correlation without delay. Although there 
was no direct correlation between cAMP and Ca 2+ , examination of cross correlations using small time 
windows showed frequent switching behavior from negative to positive and vice versa. Furthermore, 
spatially asymmetric cAMP and cGMP signals in freely deviating growth cones were visualized directly. 
These results indicate that we succeed in relating second messenger crosstalk to growth cone deviation and 
extension, and also indicate the possibility of predicting axon guidance from this second messenger 
crosstalk. 

During neural circuit formation, the growth cone navigates axons toward their targets via interactions with 
various guidance cues 1,2 . These cues activate receptors at axonal tips and change the levels of second 
messengers (Ca 2+ , cAMP, and cGMP) that play crucial roles in axon pathfinding 3 ' 4 . Imaging and optical 
tools have revealed the importance of spatially asymmetric Ca 2+ and cAMP elevation in the growth cone for the 
turning response 5 . Previous reports suggested second messenger interactions in the growth cone and neurite; 
cAMP levels affect Ca 2+ spike frequency 6 ' 7 and cAMP responses to stimuli mirror those of cGMP 8 . Moreover, the 
turning responses of the growth cone to the attractive guidance cue netrin 1 depend on the cAMP:cGMP ratio 9 . 
Although cAMP, cGMP, and Ca 2+ are key players in axon guidance, crosstalk between these second messengers in 
the growth cone has not been analyzed directly. Recent studies have reported advanced fluorescent imaging 
techniques using genetic tags and fluorescence resonance energy transfer (FRET) 10 " 12 . Therefore, to visualize 
cAMP, cGMP, and Ca 2+ dynamics in the growth cone, we performed dual FRET measurements with a single 
excitation light using Epacl -camps 13 , red cGES-DE5 n , and SapRC2.12 u , respectively. 

Results 

A confocal microscope-based dual FRET system was set up (Supplementary Fig. 1). Image processing and 
analysis, including linear unmixing, region of interest (ROI) extraction, and ROI-by-ROI FRET ratio calculations, 
were executed semiautomatically using original image analysis software written in MATLAB (Supplementary 
Figs. 2, 3). It is difficult to measure the levels of a pair of second messengers in extending thin growth cones with 
high motility; however, our imaging and post-processing system enabled the simultaneous measurement of 
cAMP and cGMP, and also cAMP and Ca 2+ in extending growth cones. We prepared dissociated cell cultures 
of dorsal root ganglion (DRG) neurons from embryonic rats (E18-19) expressing two FRET sensors, and 
performed dual FRET measurements (Supplementary Fig. 4). 

A negative correlation between cAMP and cGMP levels in non-stimulated freely extending growth cones was 
visualized directly by simultaneous imaging (Fig. la-c). The rate of axon extension was unaffected by the 
excitation laser and the expression of the FRET sensors (Supplementary Fig. 5). As expected, these results 
support previous studies using pharmacological stimuli 8,9 , and indicate that cAMP and cGMP are regulated 
reciprocally via an endogenous constitutive mechanism. Additionally, in the experiment to acquire images every 
20 s, there was no delay in their correlation in freely extending growth cones (Fig. Id). This negative correlation 
may be caused from ensembles of the phosphodiesterase family that have cAMP- and cGMP-selective character- 
istics 81415 . In contrast, there was no correlation between cAMP and Ca 2+ (Fig. le-g). However, examination of 
cross correlations using small time windows showed frequent switching behavior from negative to positive and 
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Figure 1 | Correlations between cAMP, cGMP, and Ca 2+ in growth cones, (a) Representative fluorescence image of an Epacl -camps and red cGES-DE5 
co-expressing growth cone. Lines drawn as cyan, green, yellow, and orange indicate the cell edge at 0, 20, 40, and 60 min after the onset of image 
acquisition, respectively. Scale bar represents 10 (am. (b) Spontaneous changes of AR/R average of ECFP/EYFP (cAMP) and T-Sapphire/dimer2 (cGMP) 
averaged within the ROIs shown in (a), (c) Negative correlation between cAMP and cGMP levels in the growth cones (correlation coefficient is —0.71, 
561 1 time points from 13 growth cones), (d) Cross correlation function between the spontaneous dynamics of cAMP and cGMP in the growth cones (n = 
13) showed a negative correlation without delay. Data are presented as means ± s.e.m. (e) Representative fluorescence image of an Epacl -camps and 
SapRC2.12 co-expressing growth cone. The colored lines are the same as in (a). Scale bar represents 10 (am. (f) Representative time course changes in 
spontaneous AR/R averag e values of ECFP/EYFP (cAMP) and dimer2/T-Sapphire (Ca 2+ ) averaged within the ROIs shown in (e). (g) No correlation 
between cAMP and Ca 2+ levels in the growth cones (correlation coefficient is —0.19, 8529 time points from 27 growth cones), (h) Heat map showing 
representative temporal changes in the cross correlation function between the short-term dynamics of cAMP and Ca 2+ shown in (f) (see Methods). 



vice versa (Fig. lh). The mechanism underlying this result may be 
positive feedback regulation of Ca 2+ and cAMP 6,7 . 

Spatially asymmetric cAMP and cGMP signals in freely deviating 
growth cones were visualized directly (Fig. 2a-d). There was no 
significant change in the right/left cAMPxGMP ratio (cAMP/ 
cGMP right /i e f t ) in freely extending growth cones. However, significant 
cAMP/cGMP right/left elevations were observed in the turning-side of 
the growth cones at 9.7 ± 0.5 min before growth cone angle change. 
We observed significant increases in the maximum angle change 
within 15 min after significant cAMP/cGMP right/left elevation 
(Fig. 2e). Although cytosolic cAMP/cGMP affects growth cone turn- 
ing responses to netrin l 9 , little is known about the spontaneous 
dynamics of both cAMP and cGMP. Our result suggests that endo- 
genous changes in cAMP and cGMP affect the growth cone turning 
angle. Consequently, our result supports a recent study that demon- 
strated cAMP transients driven by caged compounds steer axons 7 . 
The observed delay may reflect the signal transduction of cytoskeletal 
protein reorganization 16 . 

We next investigated the relationship between second messengers 
and axon extension. The cross correlation between cAMP and Ca 2+ 
switched frequently (Fig. lh). At the timing of switching from pause 
to extension, the cross correlation switched to negative (Fig. 3 and 
Supplementary Fig. 6). Although no common correlation was 
observed before extension, marked negative correlations were 
observed at the onset of and after extension. The rate of axon out- 
growth is reportedly inversely related to the frequency of endogenous 
Ca 2+ transients in growth cones 17 . Although we cannot discuss the 
frequency of endogenous Ca 2+ transients because of the long interval 
used for image acquisition, the dual FRET method, with a 20-s inter- 
val, revealed the relationship of cAMP and Ca 2+ to growth cone 
motility. A mathematical model based on CaMKII/CaN switching 
has been advocated and involves the upstream regulators, cAMP and 



Ca 2+1819 . The determinants of the turning response are not only Ca 2+ 
spike amplitudes but also resting Ca 2+ levels; thus, it is very import- 
ant to monitor resting Ca 2+ levels to predict the growth cone turning 
angle. 

Discussion 

During formation of specific neural connections, growth cones nav- 
igate axons toward their targets via physical and chemical interac- 
tions with environmental cues. Chemotactic substances, mechanical 
barriers, substratum adhesive substances are vital keys for neural 
circuit formation. Activation of receptors at axonal tips by these cues 
may generate signals that are transmitted back to the growth cone 
where they are transduced into changes in its motility. 

Ca 2+ signal plays a crucial role in axon pathfinding. Asymmetric 
Ca 2+ signal induces activation of CaMKII and CaN; leading to asym- 
metric cytoskeletal protein reorganization and endocytosis/exocyto- 
sis, which steers the growth cone toward attractive guidance cue 4 . 
cAMP and cGMP are known as the key regulators of cytosolic Ca 2+ 
through IP 3 R and RyR modulation. Previous reports have shown 
that the growth cone turning responses to attractive guidance cue 
netrin 1 depends on the ratio of cAMP and cGMP 9 . In this study, 
membrane permeable analogues of cAMP and cGMP were used for 
investigating the turning behavior of the axon, and direct measure- 
ment of intracellular second messengers has not been done. As such 
both cAMP and cGMP are key players in axon guidance, however 
these dynamics in the migrating growth cone are poorly known, 
especially in non-stimulated and freely extending axons. 

Our results describe the correlations between cAMP and cGMP, 
and also cAMP and Ca 2+ , between the asymmetric cAMP/cGMP 
signal and the growth cone deviation angle among cAMP, Ca 2+ , 
and axon extension. Simultaneous dual FRET imaging directly 
revealed their correlation, and is a suitable method to investigate 
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Figure 2 | Correlations between asymmetric cAMP/cGMP elevation and the deviation angle of growth cones during free extension, (a) Time-lapse 
fluorescence images of a freely extending growth cone co-expressing Epacl -camps and red cGES-DE5. The Roman numerals (i-v) indicate the time 
points shown in (b). (b) Time course of growth cone angle changes (see Methods). Right and left turns are defined as positive and negative values, 
respectively. The Roman numerals correspond those in (a). The orange horizontal line indicates the period when the maximum angle changes were 
calculated within 15 min after significant cAMP/cGMP right/left elevation (asterisk). A significant elevation was defined as a change that was 2.58-fold more 
than the derivative of cAMP/cGMP right/left and also cAMP/cGMP left/right . The cyan line denotes the latency from a significant increase to the maximum 
angle change, (c) Spontaneous changes in the right/left ratios of cAMP/cGMP, cAMP, and cGMP (see Methods), (d) Cumulative distance during the 
observation of the growth cone in Fig. 2a. (e) Freely extending growth cones deviate toward the cAMP/cGMP elevated side. The maximum growth cone 
angle changes within 15 min after significant cAMP/cGMP elevation at the turning side were higher than in randomly selected cones (80 turnings from 14 
growth cones). The randomly selected maximum growth cone angle changes were calculated as follows: select a time point from the overall time course 
and choose the maximum from within 15 min after the time point. Data are presented as means ± s.e.m. *P < 0.001. 



spatiotemporal crosstalk between two second messenger signals. 
Additionally, image processing allows us the complex analysis of 
growth cones with spatially asymmetric signal transduction cascades. 

Methods 

Cell culture. Primary DRG neurons dissected from El 8 rats (Charles river) were 
dissociated with 1 mg/mL trypsin ( Sigma- Aldrich) in phosphate-buffered saline. All 
animal procedures were approved by the Ethics Committee of Keio University 
(permit number, 09106-(1)). The dissociated cells were resuspended in Neurobasal 



medium without phenol red (Life Technologies) supplemented with 2% B27 (Life 
Technologies) and 50 ng/mL nerve growth factor (Alomone Labs). The cells were 
transfected with plasmids coding Epacl -camps and red cGES-DE5, and also Epacl - 
camps and SapRC2.12 by electroporation using Neon (Life Technologies). The cells 
were plated on a dish coated with 0.01 mg/mL poly-D-lysine (Life Technologies). The 
cells were used after a minimum of 24 h in culture. 

Dual FRET measurements. Epacl -camps is fluorescent indicator for cAMP based on 
the c AMP -binding domain of Epac and exhibits FRET that decrease on addition of 
cAMP 13 . Red cGES-DE5 is cGMP sensor based on cGMP binding domain from a 
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Figure 3 | Correlations between cAMP, Ca 2+ , and extension in freely extending growth cones, (a) Time-lapse fluorescence images of a freely extending 
growth cone co-expressing Epacl -camps and SapRC2.12 at the time points shown in (b). The Roman numerals (i-vi) indicate the time points shown in 

(b) . Scale bars represent 5 um. (b) Cumulative distance during observation of the growth cone in Fig. 3a. The Roman numerals correspond to those in (a). 

(c) Time course of spontaneous changes in the AR/R average values of ECFP/EYFP (cAMP) and dimer2/T-Sapphire (Ca 2+ ). The colored horizontal bars 
indicate the 10-min time windows used to calculate the correlations shown in (d). The blue, cyan, and other colored bars denote before extension, start of 
extension, and after extension, respectively. Cross correlations between cAMP and Ca 2+ become negative at the change from pause to extension. 

(d) Temporal changes in the cross correlation function between the short-term dynamics of cAMP and Ca 2+ before and after extension. The colors of the 
lines indicate the timing of the cross correlation calculations denoted in (c). (e) Heat map showing the averaged sequential time course of the cross 
correlation function between the short-term dynamics of cAMP and Ca 2+ . The growth cones started to extend at 0 min (15 extensions from 11 growth 
cones). The start of extension was defined as when the growth cones extended 5 um in 10 min. 



phosphodiesterase 5 and also exhibits FRET as Epacl -camps 11 . SapRC2.12 is Ca 2+ 
indicator constructed by fusion of T-Sapphire, CaM, M13 fragment from myosin 
light chain kinase, and dimer2. SapRC2.12 exhibits FRET that increase upon increase 
inCa 2+u . 

Dual FRET imaging experiments were performed using a laser scanning con- 
focal microscope system (Olympus) mounted on an inverted microscope 
(Olympus) with a X 100 oil immersion objective. Both donors of the FRET 
sensors were simultaneously excited using a diode laser at 405 nm, and the 
excitation was split using a beam splitter 20/80, the emission was split using 3 
dichroic mirrors (490, 510, and 560 nm), and the emissions were acquired using 2 
grating-based detectors (460-490 nm and 490-510 nm) and 2 filter based detec- 
tors (535-565 nm and 560-660 nm). Fluorescence images were acquired every 



20 s. Neurons were maintained at 37°C in a humidified atmosphere of 5% C0 2 
(Tokai Hit). 

Image processing. Images were analyzed using original image processing and 
analysis software written in MATLAB (Math Works). Acquired images were 
smoothed using a 3 X 3 median filter to reduce noise and the background was 
subtracted. To obtain the fluorescence signal of each region, ROIs were calculated as 
described in Supplementary Fig. 2. To obtain the fluorescence signal of each 
fluorescent protein, an unmixing processing was applied to the fluorescence signals 
averaged within the ROIs 11 . Finally, the FRET signals were calculated from the 
unmixed fluorescence ROI-by-ROI. The growth cone angle and movement were also 
calculated as described in Supplementary Fig. 3. 
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Statistical analysis. Cross correlations and all statistical analyses were performed 
using MATLAB (Math Works). To investigate the correlation between cAMP and 
Ca 2+ , we calculated the temporal changes in cross correlation between the short-term 
dynamics of cAMP and Ca 2+ . cAMP and Ca 2+ levels in 10-min time windows 
normalized to their average were used as input for the cross correlation calculations. 
Significant differences were determined using Student's f-test. *P < 0.001. 
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